Number of spanning clusters at the high-dimensional percolation thresholds 
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A scaling theory is used to derive the dependence of the average number (fc) of spanning clusters at 
threshold on the lattice size L. This number should become independent of L for dimensions d < 6, 
and vary as log L at d = 6. The predictions for d > 6 depend on the boundary conditions, and the 
results there may vary between and L". While simulations in six dimensions are consistent 

with this prediction (after including corrections of order log(logI/)), in five dimensions the average 
number of spanning clusters still increases as log L even up to L = 201. However, the histogram P{k) 
of the spanning cluster multiplicity does scale as a function of kX{L), with X{L) = 1 + const/L, 
indicating that for sufficiently large L the average (k) will approach a finite value: a fit of the 5D 
multiplicity data with a constant plus a simple linear correction to scaling reproduces the data very 
well. Numerical simulations for d > 6 and for d = 4 are also presented. 
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I. INTRODUCTION AND THEORY 

We are interested in site percolation on a finite hy- 
percubic lattice in d dimensions, of linear size L. The 
number of spanning clusters when L is of the order of 
the percolation correlation length, ^, away from the per- 
colation threshold Pc, has been discussed originally in 
. The purpose of this paper is to discuss the average 
number of spanning clusters, (fc), at p = pc, when ^ is in- 
finite, and all the critical quantities behave as powers of 
L. As we discuss below, the number of spanning clusters 
has been the topic of much discussion in the literature. 
For example, although intuitively one would think that 
in two dimensions there exists only one spanning cluster, 
in fact (for appropriate boundary conditions) there is a 
whole distribution of the sizes of such clusters 0|. In 
addition, we show below that for d < 6 the average num- 
ber of such clusters is finite, in direct relation with the 
validity of the hyper-scaling relations among critical ex- 
ponents. For d > 6 hyper-scaling is violated, dangerous 
irrelevant variables must be introduced, and the result 
for (k) becomes ambiguous, depending on the boundary 
conditions. In fact, the theory for that case is yet incom- 
plete, leaving an open challenge for future research. 

We start with a theoretical discussion. The percola- 
tion "order parameter" Pspan is usually defined as the 
probability that a site belongs to any spanning cluster 



(i.e. to the union of all spanning clusters). For a finite 
sample, Pspan is related to the cluster size distribution 
function ns{p, L) (defined as the average number per site 
of clusters containing s sites) via the sum rule 3] 

Pspan{p,L) ^p~^sns{p,L), (1) 

s 

where the sum over s excludes all clusters which span 
the lattice. The exact details depend on the definition of 
"spanning" , e.g. along how many directions should the 
cluster connect opposite faces of the hypercube. How- 
ever, these details do not matter for the scaling argu- 
ments presented below. 

For a finite lattice, the sum in Eq. ^ goes up to 
Smax{L), which is of the same order as the average mass 
of a single spanning cluster, s{L). Since the total mass 
of all spanning clusters is given by L''' Pspan , this implies 
that the average number (k) of spanning clusters is given 

by 

(fc) « Pspan /s{L). (2) 

This relation should hold at all dimensions. The propor- 
tionality constant in Eq. ^ (which results, among other 
things, from dividing averages rather than averaging the 
ratio and on the detailed definition of "spanning" in a 
finite sample) may depend on the boundary conditions 
and on other details. 
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One way to derive s{L) is to use the pair connectivity 
function, G(r). This function yields the probabihty that 
a site at distance r from the origin is connected to the 
same cluster as the origin. Assuming that the site at the 
origin belongs to any spanning cluster, with probability 
Pspan, the density profile of the cluster is given by p{r) = 



G{r)/Ps, 



Lensity 

m 



. Thus, 



S{L) CX / d\ G{r)/P,pan = S{L)/ Pspan, (3) 

Jo 

where S{L) — cfir G{r) is proportional to the mean 
cluster size. 



S CX ^S^Usip, L)/^; 



(These two quantities are equal in the thermodynamic 
limit, L — > oo.) Thus, we come to the fundamental rela- 
tion 



(fc) « L^Pl^JSiL). 



(4) 



A naive finite size scaling theory would predict that 
at p = one has P^pan oc P-l^/" and S{L) oc L^l" 0. 
Thus, one concludes that 



(5) 



For d < 6, hyper-scaling implies that dv — 2/3-1-7, ^"^^ 
therefore (fc) is asymptotically independent oi L. In fact, 
the combination of amplitudes which appears in Eq. (0) 
is universal, depending only on the type of boundary con- 
ditions • As we discuss below, numerical estimates for 
{k) indeed approach a constant for d < 4. However, data 
for d = 5 require further discussion. Although not a ma- 
jor purpose of this paper, our discussion below should 
also serve as a warning and as a guideline for future nu- 
merical simulations in such high dimensions. In fact, 
some data remain ambiguous unless one uses available 
theoretical information, or unless one performs simula- 
tions on much larger scales than presently possible. 

At d = 6, many power laws are modified by logarithmic 
corrections. In fact 0, 



5(L) ocL2(lnL)i/2i. 
Using our basic result, Eq. we thus find 



r6p2 



(X In L, d = 6. 



(6) 



(7) 



Further analysis shows that the coefficient of proportion- 
ality here (and also the result for (k) for d < 6) is a uni- 
versal number 0. Corrections to the above result will 
involve ln(lnL-|- const). Fig. shows that the data in 
six dimensions, which were fitted by the square of log L 
in 0, can also be fitted to a simple logarithm plus a log 



log L correction to scaling, (fc) = A In L -I- i? In In L + . . . . 
The two data sets refer respectively to free boundary con- 
ditions (FBC) and to mixed boundaries (MBC), i.e. he- 
lical in d — 2 directions and free in the remaining two. 
In both cases we obtain a good fit for L > 10, so that 
the trend appears independent of boundary conditions, 
as one would expect. The coefficients of the fits are 
A = 11.9, B = -17 (MBC) and A = 2.38, B = -4.33 
(FBC). 




FIG. 1: Average number (k) of simultaneously spanning clus- 
ters at the six-dimensional percolation threshold, for FBC 
(*) and MBC (+) (see text). The lines are fits to (k) = 
A\n{L) + _Bln(lnL); we used the data of Q plus additional 
simulations for small L with the number of samples variable 
from 1000 to 50000. 

For d > 6, hyper-scaling is broken, and one has 
the mean-field exponents /3 = 7 = 2zy = l. A sim- 
ple scaling in which ^ is replaced by L would give 

Pspan L 



S{L) CX L^, s{L) CX and 
(fc) CX L'^- 



(8) 



This simple result, which already appeared in 0, has 
also been proved by Aizenman [2| for the case of bulk 
boundary conditions, where a spanning cluster connects 
two opposite faces of the box of size L under the condi- 
tion that sites in the box can also be connected by paths 
outside the box. 

However, other forms of scaling may be possible, which 
may depend on the boundary conditions. We will discuss 
the different scaling approaches within the framework of 
the renormafization group (RG) theory (e.g. Ref. .7]). 
This theory is conveniently discussed using the "free en- 
ergy" F, equal to the generating function of the cluster 
distribution function, 



F{p, h, L) = ^ns{jp, L)e 



— sh 



(9) 



which is related to the free energy of the g-state Potts 
model for q ^ 1. Pspan and S{L) are then the first and 
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second derivatives of the singular part of F with respect 
to h. After £ RG iterations, this singular part becomes 



f{p,h,w,L) 



f{t{£),h{e),w{£),L/e' 



(10) 



where t = p — Pc- Here, w represents the probability for 
having a three-fold vertex at a site on a cluster. Although 
irrelevant in the RG sense, this variable must be included 
in the analysis, since the "free energy" may depend on 
it in a singular way, causing the breakdown of hyper- 
scaling. This is why w is called "a dangerous irrelevant 
variable" . For a finite sample close to pc, L iteration 
until — L yields the L-dependence of the quantities 
discussed above. 

For d < 6, w{£) approaches a finite fixed point value, 
and one ends up with hyper-scaling and with the conclu- 
sion that k approaches a finite constant. For d > 6, one 
ends up with 

fit, h, w, L) = L-'^f{tL^, hL'^/^+\wL^-'^/^, 1). (11) 

However, as stated above, now w turns out to be a "dan- 
gerous irrelevant variable". One way to see this is to 
consider the infinite sample. After eliminating the fluc- 
tuations, / is given by the minimum of the Landau free 
energy, 

wP^-hP], (12) 



f = mm\tP^ 
p 



while the order parameter Pspan is equal to the value 
of P which minimizes this free energy. Replacing P by 
Qw^^+^Z^, with an arbitrary exponent x, it is easy to see 
that / obeys the scaling relation 

fit, h, w) = w;-2+^/(tu;-^/3^ hw^-^'^/^, 1). (13) 

In particular, the minimization with respect to Q now 
yields the equation 3Q^ + 2Qtw-^^^ - hw^-^""/^ = 0, 
leading to the scaling form 



fit,h,w) = w-'t^gihw/t'). 



(14) 



which is completely independent of the arbitrary expo- 
nent X. This ambiguity stems from the fact that the 
Landau free energy is calculated for an infinite system. 

Unlike the above result for the infinite system, the ex- 
ponent X does persist for finite samples. In fact, the 
theoretical predictions for (fc) depend crucially on x, and 
this remains an open challenge for future research. In 
this case, combining the scaling with L from Eq. Hll() 
with Eq. H13|l . we end up with 

fit, h, W, L) = y;^-2i-6+3x-d:r/2 ^ 

and the value of x must follow from the boundary con- 
dition, which breaks the scale invariance reflected in Eq. 

Taking derivatives with respect to h, and setting 
t = {), this implies that 

^ span ^ : 
S OC 

siL) OC L4-2x+d^/3^ (jg) 



The naive scaling result of Eq. 
simple choice x — Q. 



is obtained with the 
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FIG. 2: Efltective exponents for the scaling with L of the span- 
ning cluster multiplicity (from f^), as a function of the spatial 
dimension d. The dashed line represents the theoretical pre- 
diction from Eq. Q. 

Fig. [2 compares Eq. ((SJ with the numerically es- 
timated effective exponents. Clearly, there exist some 
discrepancies between theory and simulations. However, 
these discrepancies do not worry us much since the lat- 
tices for d > 6 are very small in their linear dimension 
L. This is confirmed by the fact that, for c? > 6, the 
results for the exponents strongly depend on the type of 
boundary conditions one chooses. Also, our attempts to 
add corrections to scaling to the leading term of the fit 
ansatz did not improve the situation. At present, the 
numerical results cannot clearly confirm that a; = 0, as 
would be required if Aizenman's result also applies 
for our boundary conditions (note that his proof works 
under somewhat different conditions). Thus, the value of 
X remains to be determined in the future. 

A second possibly relevant case is x = 2. In this case, 
/ has the form 

fit, h, w, L) = L-^Tiw-^lhL^/^ hw-^'^'L^^'^ 1, 1), 

(17) 

yielding D = 2(i/3 and Pspan and S given by 



P> 



span 



^L-^/'^f^iw-^'W^), 



(18) 



S^L'"^f2iw-^/HL''/''). (19) 
In the limits t {), L ^ oo, these equations behave as 



Ps, 



Pspan ^ S 



(X S (X i cx ^ 



L = oo. 



(20) 



OC L 



-d/3 



(k) OC const, t = 0, (21) 
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and 



P. 



(k) oc L'^-^y, 



tL^ = const < 0, y < d/3. 



(22) 



(Remember: these equations only apply for d > 6). The 
scaling behavior of Eqs. H20l) - (|22|l for j/ = 2 is analogous 
to that found by Chen and Dohm 9] for the (j>* Ising 
model with periodic boundary conditions, and has re- 
cently been used to find an upper bound for the number 
of spanning clusters |^] . As seen from Fig. [S] our data for 
t = do not seem to agree with the prediction that (k) 
approaches a finite constant. However, our samples are 
small, and the situation for d > 6 requires more studies. 
Very recent simulations show that Pspan differs apprecia- 
bly from the relative size of the largest cluster [T^l , but 
that is not enough to conclude that the spanning cluster 
multiplicity diverges for c? > 6. Thus, the true behav- 
ior of (k) for d > 6 remains an open issue. We hope 
that the present paper will stimulate both numerical and 
theoretical discussions of this question. 



II. RESULTS OF THE SIMULATIONS 



From our Eq.Q we expect that for L large this ratio 
converges to a constant, and both our figures confirm 
this expectation. We stress that the S that we calculate 
through our simulations differs from the standard defini- 
tion of S" = J2s s'^^sip, L)/ sus by the absence of the 
denominator. The latter is smaller than f, as it is the 
density of occupied sites which belong to finite clusters, 
so the real value of the ratio L'^Pspa,i/{S (k)) would be 
smaller than the one we show in Figs. 13 and ^ On the 
other hand the functional dependence on L of the ratio 
is the same in both cases. 




We now review in detail our numerical simulations. 
The first numerical studies on this topic were performed 
by De Arcangelis 
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FIG. 3; Numerical verification of our relation Eq. JIJ, for 5D 
percolation with MBC. 

An important issue concerns the validity of our Eq. 
which is a general result that establishes a link be- 
tween the order parameter Pspan, the mean cluster size 
S and the spanning cluster multiplicity (k). Figs. El and 
2] show numerical tests of this relation. We analyzed one 
case below the upper critical dimension dc = 6, i. e. 5D 
with MBC, and one case above dc, i. e. 7D with FBC. In 
both cases we have calculated the ratio L'^Pgp^^/{S (fc)). 



FIG. 4: As Fig. El but for 7D percolation with FBC. 

Let us now concentrate on the trend of the multiplicity 
when the linear dimension L of the lattice varies. In less 
than five dimensions, Ref. found an asymptotically 
constant number of spanning clusters for both theory and 
simulation. Thus the real problem is five dimensions, 
which we discuss now. 

Fig. shows new data for sites, extending up to 
L ~ 201, the largest five-dimensional system known to 
us from direct simulations (176^ was simulated in Ising 
models 112! )• They are still compatible with the number 
of spanning clusters increasing as log L. These data use 
MBC, but similar proportionalities to log L were found 
also with free and periodic boundary conditions (PBC, 
which means helical boundaries in d— \ directions), which 
are illustrated in Fig. By looking at the data corre- 
sponding to periodic boundaries, however, it seems that 
the curve smoothly bends towards the end, which might 
indicate the beginning of a crossover. 

For comparison. Fig. shows four-dimensional sim- 
ulations with MBC. There we see for small L < 12 a 
logarithmic increase, followed by a crossover region ex- 
tending over one decade in L, and ending finally in the 
theoretically predicted plateau near L — 200. Further 
simulations with PBC (not shown) confirm this trend. 
These four-dimensional results are a nice example for the 
need of large lattices in simulations: only above L = 100 
the theory is confirmed. Comparing Figs. and |3 we 
may hope that also in five dimensions a crossover would 
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FIG. 5: Five-dimensional problem, MBC. The right-most er- 
ror bar for L = 201 comes from 10 samples, the two smaller 
error bars shown from 100 samples. Most of the data use 
1000 samples with error bars too small to be shown. There 
is no evidence for a plateau; nevertheless, the trend of the 
data is beautifully reproduced just by adding a simple non- 
logarithmic scaling correction (dashed line in the plot). 



be seen towards a plateau, if we could simulate larger 
lattices than the present world records [T3l |. 
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FIG. 6: Five-dimensional data of the spanning cluster mul- 
tiplicity at threshold for FBC and PBC (*). One still 
finds an increase with L, although for PBC it seems that the 
data curve starts bending at high L. The number of samples 
goes from 3000 to 50000 for FBC and from 10000 to 50000 
for PBC. 

There is a way, however, to clarify the issue. Instead 
of looking merely at the average number (fc) of spanning 
clusters (as in all of the above discussion), we can analyze 
the histogram P{k) of that number. Its tail for large 
k was already shown to be consistent with theory |l4| : 
- In P{k ^ oo) oc for 2 < d < 7. If (fc) would 

be size-independent, on a lattice of linear dimension L 
we would also expect P{k,L) = P{k), i.e. the histogram 



is as well size-independent. Finite-size scaling for large L 
would nevertheless allow correction factors like X{L) = 
1 + const/L, so that the histogram is indeed a function 
P{kX{L)). On the other hand, P{kX{L)) is normalized 
to some L-independent number C (we chose here C = 
1000), and therefore 



£p(fcX(L)) = C. 



(23) 



fe=0 



If we approximate the sum with the integral over fc, we 
can perform the change of variable j = kX{L), so to get 



1 



XiL) 



Y.P{])^C. 



(24) 



Equation (|24|l cannot be right as it stands, because the 
left hand side is a function of L, while the right hand 
side is not. That means that the histogram P(fc, L) is 
not a scaling function of the variable j = kX{L)^ but 
that P{k,L) = X{L)P'{kX{L)), where P' is now a 
scaling function of j. Let us check what happens to the 
average multiplicity (fc) if we assume that P(k,L) has 
the above-derived form: 



{k) = 



1 



Y.Z,kX{L)P'{kX{L)) 

Yr=,X{L)P'{kX{L)) CX{L) 



where we again made the substitution j = kX{L) in the 
sum over fc. The sum in (|25|l is independent of i, and 
we finally obtain: 



(fc)cx 



1 



x{Ly 



(26) 



so that the whole L dependence of the average multiplic- 
ity is contained in the correction factor X{L). 

In this way, we have now the chance to make a cross 
check on our data. If (fc) is indeed size-independent 
for d < 6, we should be able to find a simple non- 
logarithmic correction factor X{L) ~ 1 + const/i, such 
that (fc) cx l/X{L) and consistently P{k,L)/X{L) is 
a scaling function of the variable j = kX{L). We re- 
mark that only eventual discrepancies of the rescaled his- 
tograms at large fc can lead to infinitely many spanning 
clusters in the limit L — > oo: if, for all L, there were a fi- 
nite fc* beyond which the histograms collapse, (fc) would 
necessarily approach a constant for large L, modulo finite 
size corrections. 

We tried then to fit the multiplicity curves in four and 
five dimensions with the simple two-parameters ansatz 
aL/{L + b), to see whether we could reproduce the 
observed trends. As a first trial we took the four- 
dimensional data of Fig. [7| and fitted several portions of 
the increasing part of the curve, before the plateau. All 
our fits are quite good; besides, starting from the range 
[2 : 60], the fit is stable, i.e. we obtain the same values 
for the parameters a and b as for the fit on the full curve, 
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FIG. 7: Four-dimensional solution, MBC. The spanning clus- 
ter multiplicity increases in a wide range of L, but attains 
finally a plateau for L > 100. The dashed curve is the fit 
with the simple correction to scaling ansatz aL/(L + b). For 
each L we took mostly 1000 samples. 

within errors. This is quite interesting, because it allows 
us to predict quite precisely where saturation takes place, 
even if one analyzes values of L which lie well below the 
beginning of the plateau. The best fit curve, for which 
a = 0.78 and b = 4.2, is plotted in Fig. □ 
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FIG. 8: Four-dimensional histogram, MBC. Occurrence 
P{k, L) of samples with k spanning clusters each, versus fe, 
with the finite-size correction factor 1 -I- 4.2/L derived from 
the fit of the average multiplicity. The data sets have a vari- 
able number of samples, mostly 50000, we normalized them 
all to 1000 iterations. The nice scaling is consistent with the 
size-independence of the average multiplicity. The lattices we 
have taken are: 20*, 30"*, 40", 50", 70^*, 85^. 

As a nice confirmation of this result, we show in Fig. |S| 
the corresponding histograms P{k,L), where we use the 
value b = 4.2 derived above for the scaling correction. 
We chose again on purpose only values of L before the 
plateau of the average multiplicity. The figure gives a 



nice data collapse, as we expected. 
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FIG. 9; Same as Fig. |S| for five-dimensional data with MBC. 
One obtains again a very good scaling of the spanning cluster 
multiplicity distributions for A; > 0, by introducing the simple 
correction 1-1-17.5/L. The lattices are: 30^, 40^, 50^ 60^, 70^, 
80^ 91^ The number of samples goes from 1000 to 50000. 

We repeated the analysis for the five-dimensional data, 
starting from the puzzling curve of Fig. [3 Here we find 
that the fits are very good and stable starting from the 
very beginning of the curve: the fit parameters a and b 
are basically fixed already in the range [2 : 20]; the best 
fit (dashed line in the figure) was obtained including all 
datapoints with significant statistics, i.e. for 2 < L < 91; 
we obtain a = 3.12(6), b — 17.5(7). As one can see in 
Fig. our simple ansatz (dashed line) describes very 
well the observed behavior of the data. We notice that 
on the logarithmic scale for L our simple scaling curve 
has an inflection point, exactly as the data. This inflec- 
tion can be interpreted as a signal of a possible crossover 
from an initial logarithmic increase of the multiplicity to 
a successive convergence to a plateau; we now see that it 
is instead a natural feature of our scaling ansatz. Notice 
that the correction term is more important than in four 
dimensions (17.5 vs 4.2). That means that the data con- 
verge much more slowly to the plateau, and explains why 
we could not see a saturation even at L = 201 (although 
the argument can be reversed). Indeed, for a given L, 
the ratio r of the multiplicity to the asymptotic plateau 
is r = L/{L + b). In 4D, r = 96% for L = 100 and 
r = 98% for L = 200; in 5D one would obtain 85% and 
92%, respectively. In order to "see" the plateau as we 
do in four dimensions, we would need to go to L ~ 800! 
To check the consistency of the picture in 5D we studied 
the scaling of the histograms P{k, L), with the correction 
constant b = 17.5 that we determined above. The result 
is illustrated in Fig. the scahng is quite good, except 
eventually for k — (but scaling laws seldom hold for 
small integers) and at the very end of the tail, where the 
statistics is too low and there are relevant fluctuations of 
the data points. 
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of L (from 10 to 80); all our fitting curves include a cor- 
rection term, for which we fixed the value of the exponent 
uj to the estimate 0.53 given in Taking into account 
this correction, the finite size scaling fits are remarkable 
for all percolation variables if we use the PBC data, for 
which we get = 1.45(2) and 7/1/ = 2.08(2). For 
MBC the fits are also very good, but a bit worse as far 
as the and the values of the exponents are concerned 
(/3/t/ = 1.45(2), j/i^ = 2.10(2)); for the FBC the fits 
are not so good and the values of the exponents not in 
agreement with expectations, though quite close. Also 
series expansions [l6j , which are independent of any lat- 
tice size, gave in five dimensions the usual exponents in 
agreement with expectations, without indications of par- 
ticular difficulties. 



FIG. 10: Same as Fig. ^ but with FBC. The correction is 
1 + 4/L, the lattices are: 20^ 30^ 40^ 50^ 60^ 70^ The 
number of samples goes from 3000 to 50000. 
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III. CONCLUSIONS 

In summary, we have derived the scaling behaviour at 
threshold of the average number (k) of spanning clusters 
with the linear dimension L of the lattice, for any space 
dimension d. Below the upper critical dimension dc = 6, 
(k) should approach a constant when L — s- 00, for d = 6 
it should increase as log L, and for d > 6 it could increase 
as L''^^, but could also approach a constant (depending 
on boundary conditions and yet unknown theoretical de- 
tails). While the latter conclusions might seem just a 
confirmation of previous results on the topic, our work 
highlights two new important issues: 

• the possibility of other scaling behaviors of (k) 
above the upper critical dimension, which possibly 
depend on the boundary conditions; 



FIG. 11: Same as Fig. ^ but with PBC. The scaling is re- 
markable, probably due to the better statistics; the necessary 
correction is 1 + 20/ L. The lattices are: 30^ 40^ 50^ 60^ 
70^ 80^ The number of samples goes from 8000 to 50000. 



Finally, we analyzed the two other data sets in five 
dimensions, i.e. the ones relative to FBC and PBC. In 
both cases we found that our picture works: we could 
find a scaling correction X{L) = 1-1- const/L such that 
both the average multiplicity and the histograms show a 
clean scaling. The histograms are shown in Figs. 1101 and 
111! the scaling is again very good for fc > 0, with some 
fluctuations at the end of the tail which are likely due to 
the low statistics of those points. 

We also checked for five dimensions whether other per- 
colation quantities behave unusually, and found that they 
do not. The size of the largest cluster at the percola- 
tion threshold varies asymptotically as L'^~f^/'^ and the 
"mean" cluster size as L^/'^, where j3/v ~ 1.46, 7/:/ ~ 
2.07 are expected !3| in five dimensions. We found that 
corrections to scaling play an important role in this range 



• the relevance of corrections to scaling, which may 
affect the scaling behaviors up to two-three orders 
of magnitude in L. 

Our numerical investigations confirm that the multi- 
plicity indeed converges to a constant for c? < 6. For 
the case d = b, where we do not see a plateau even for 
the largest L we have taken, a simple linear correction to 
scaling is able to reproduce the observed data pattern. In 
six dimensions the results can be made consistent with 
theory by adding a logarithmic finite-size correction; in 
more than six dimensions both the data on the multi- 
plicity k and those on the order parameter Pgpan and the 
mean cluster size S lead to very different values of the 
finite-size scaling exponents for different sets of boundary 
conditions, and we are not able to derive reliable conclu- 
sions. So, if our numerical evidence below dc is con- 
clusive, to close the issue above dc simulations at much 
larger L seem to be necessary. 
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